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ABSTRACT 

We study the clustering properties of the first galaxies formed in the Universe. We 
find that, due to chemical enrichment of the inter-stellar medium by isolated Population 
III stars formed in mini-halos at redshift z > 30, the (chronologically) first galaxies are 
composed of metal-poor Population II stars and are highly clustered on small scales. 
In contrast, chemically pristine galaxies in halos with mass M ~ 1O 8 M0 may form at 
z < 20 in relatively underdense regions of the Universe. This occurs once self-enrichment 
by Population III in mini-halos is quenched by the build-up of an H2 photo-dissociating 
radiative background in the Lyman- Werner bands. We find that these chemically pris- 
tine galaxies are spatially uncorrelated. Thus, we expect that deep fields with the James 
Webb Space Telescope may detect clusters of chemically enriched galaxies but individ- 
ual chemically pristine objects. We predict that metal-free galaxies at 10 < z < 15 have 
surface densities of about 80 per square arcmin and per unit redshift but most of them 
will be too faint even for JWST. However, the predicted density makes these objects 
interesting targets for searches behind lensing clusters. 

Subject headings: cosmology: theory — early universe — galaxies: high-redshift 



1. Introduction 



The concept of first galaxies is not as easy to define as that of first stars and it's partly a 
matter of convention. A first, somewhat arbitrary, decision, is the mass scale where an object 
becomes a galaxy. A common convention in the context of the first galaxies is describing galaxies 
as those objects forming within dark halos massive enough that their gaseous content can cool 
by atomic hydrogen (Lyq cooling), i.e. characterized by a virial temperature greater than 10 4 K 
(jrlicotti et alJbood : I Wise & Abelbood : breif et alJbood ). The advantage of this definition is that 
it is based on a physical process and makes galaxy formation relatively independent of the presence 
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of a Lyman- Werner (LW) background that could hamp er the formation o f Population III stars 



in mini-halos by photo-dissociating molecular hydrogen (ILepp fc Shul] 



Haiman et al.lll997l ; iMachacek et al.ll2003l ; iTrenti k Stiavellil 12009; 



cp 

li : 



1983 



Tee mark et all Il997l : 



Trenti et al 



2009^ . 



Once we have defined "galaxy" we need to define "first" . The usual assumption is that "first" 
refers to a chronological sequence of events, so that the first galaxies are those that formed at the 
earliest redshifts. It is sometimes assumed that the chronologically "first" galaxies will also be the 
least chemically evolved ones. In fact, these objects are most likely to form in highly clustered 
areas and they are also ver y likely to have been polluted by m etals produced by Population III 
stars forming in mini-halos (jWise &: Abelll2008l ; iGreif et al.ll2008l ). Thus, these chronologically first 
galaxies would not have a primordial chemical composition. However, we can imagine another class 
of galaxies that would form after a LW background is in place and would then not be polluted 



by stars in mini-halos (jjimenez &: Haiman 



2006 



Trenti et al.1 |200J; l.lohnson et al.l l2008h . These 



objects could conceivably form stars out of pristine primordial gas and would then be galaxies 
made of Population III stars. Assuming that they can form Population III stars in significant 
numbers, these objects would be chemically less evolved (i.e. "younger") despite forming later in a 
chronological sense. 

This paper is devoted to exploring the validity of the two scenarios described above and to 
studying the redshift distributions and clustering properties of these two classes of objects. These 
quantities are relevant for planning surveys to study the first galaxies with the James Webb Space 
Telescope and other (future) facilities, such as thirty-meter class telescopes used for imaging behind 
lensing clusters. 

The structure of the paper is the following. In Sec. 2 we discuss our model for galaxy formation 
during the Dark Ages, in Sec. 3 we present our results on star formation rates and clustering. Sec. 
4 concludes with prospects for future detections and speculates on the poss i bility that metal-free 
galaxies have already been observed behind a lensing cluster by lStark et al.l ((2007). 



2. Galaxy formation during the reionization era 



Our met hod is based on a combin a tion of nume r ical s i mulations and anal y tical c onsiderations 



develo ped b vlTrenti fc Stiavellil ( 2007TI . ITrenti et al.1 (|2008T ). ITrenti Stiavellil (|2009|), ITrenti et al 



(2009), and ITrenti Shulll (|2010l ). Here we briefly summarize the main properties of our model. 



We determine the dark matter h alo formation rate and m erging history usin g either extended 
Press-Schechter modeling based on the lSheth &: Tormenl (119991) mass func tion (see lTrenti &: Stiavelli 
2009) or high-resolution cosmological simulations (see lTrenti et al.ll2009l ). We populate metal-free 
halos with stars by taking into account the cooling time scales for atomic and mo lecular h ydrogen. 



Stiavelli 



inclu ding also radiative feedback due to LW radiation (jTrenti Stiavellil 120091 . see also 
2008). The minimum halo mass required for metal- free star formation from our model is shown as 
solid black line in left panel of Fig. [TJ Metal-free star formation transitions from mini-halos with 
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virial temperature T v i r ~ 10 3 K, where molecular hydrogen cooling is efficient at z > 30, to more 
massive halos (T v i r > 10 4 K) at z < 15, where Ly-a cooling is possible independently of the radi- 
ation background. The minimum mass we require for star formation is consistent with the results 
from detailed hy drodynamic simulations o f metal free star formati on in the presence o f a radiative 
LW background (jO'Shea Sz Normanll2008r). as discussed in detail inlTrenti et al.l (120091). If a ha lo is 



able to cool and if T vir ~ 10 3 K (jAbel et al 



2002 



Yoshida et al. 



2003 



O'Shea fc Normanll2007h . we 



assume it forms a single, massive Pop III star drawn from a Salpeter initial mass function (IMF) 
in the m ass range 50-300 M^. When appropriate the star explodes as a Pair Instability Supernova 
(PISN) (jHeger fc Wooslevl [2002J; IScannapieco et al.l 120031 ) . A constant fraction of metal- free gas 
(/* « 5 x 10~ 3 ) is converted into stars when Lya cooling is possible in halos with T v \ v > 10 4 K at 
lower red shift. In this case, w e assume a lower characteristic mass (O(30)Mq) for Population III 
stars (see lYoshida et al.l l2006). 



The transitio n from metal-free to metal enriched star for mation happens at a critical metallicity 



10~ 3 - 5 Z Q bromm fc Larsonlhooi l ISmith et al.1120091 ) and might be as low as Z CT it ~ 10 6 Z Q 



in the presence of dust ([Schneider et al.ll2006l ). A single PISN provides enough metals to trigger 
the transition to metal-enriched star formation in halos with 10 4 K (that is with mass 

M < 3 X 10 7 at z > 10 - see left panel of Fig. [Q). In fact, several studies sh ow that these halos 



may even reach a metallicity well abo ve Z cr n as a result of PISN enrichment (jWise &; Abel 



Karlsson et al. 



2010 



Greif et al J 12010 ). 



2008 



For metal enriched halos, we assume that a fraction /, ~ 5 x 10~ 3 of gas is converted into 
stars with a Salpeter IMF in the mass range 1-100 Mq. In the analytical model, metal enrichment 
is based on a statistical approach making use of the pro bability that progenitors of a halo will have 
formed a star before the halo collapses, as described in lTrenti &: Stiavellil (|2007l ). In the numerical 
simulations, we regularly save full snapshots of the system and use them to construct a detailed halo 
tree to determine whether a halo is self-enriched. Through simply counting photons we track the 
establishment of an average radiative background, i.e. our models do not include detailed radiative 
transfer. 

The numerical simulations have been run using the Particle Mesh Tree Code Gadget-2 (Springel 
2005). We adopt the fifth year WMAP concordance cosmology with = 0.72, Vt m = 0.28, 
Qf, = 0.0462, erg = 0.817, n s = 0.96, h = 0.7. We start our main simulation at redshift z=199 using 
a box with edge 7/i _1 Mpc, N=1024 3 dark matter particles, a mass resolution of 3.4 x 10 M and 
a force resolution of 0.16 /i -1 kpc. Halos are identified with a friend-of-friend halo finder (Davis et 
al. 1985) using a linking length equal to 0.2 the mean particle separation. In addition to in situ 
enrichment by progenitor mini-halos, in the numerical simulations we also consider enrichment by 
winds which is not included in the analytical model. We assume an outflow s peed of 60 km s -1 
from halos that contain metal-enriched galaxies, that is with T v i r > 10 4 K (see iTrenti et al.l 12009 
for further details). The resulting wind enrichment is consistent with wh at is obtained using a 



non-cosmological Sedov- Taylor model (e.g. Eq. 8 in iTumlinson et al.l 120041 ) . namely bubble sizes 



of < 150 h x kpc at z > 6. Gas polluted by winds in dark matter halos that have not been self- 



- 4 - 



enriched is also likely to reach a metallicity Z ~ 10~ 3 ' 5 Z Q sufficient to trigger the transition to 
extremely metal poor (EMP) star formation. In fact, the typical distance traveled by winds before 
they encounter a non-self enriched halo is of the order of ~ 50/i _1 kpc, sufficient to reach the critical 
level of me tal pollution for m etal outflows originating from a dwarf galaxy (see Section 2.1 and 
Figure 3 in lTrenti et alll2009l ). 



Results 



3.1. Galaxy Formation Rate 



In Fig. [2] we show, as a function of redshift, the rate of proto-galaxies forming per year and 
made of chemically enriched stars (marked as PopII, dashed red line) compared to that of proto- 
galaxies made of Population III stars (solid blue line). The figure is derived from the analytical 
model for galaxies with virial temperature between 10 4 and 2 x 10 4 K - roughly corresponding to 
< 10 8 Mq (see left panel of Fig. Q]) - and confirms our expectation that chemically enriched galaxies 
form at earlier redshifts than galaxies made of population III stars. The analytical model does 
not include the effect of enrich ment by winds wh ich will lead to ~ 30% lower numbers of pristine 
objects at z ~ 10 (see Fig. 4 in lTrenti et al.l 12009 ) . 



In our model, chemically enriched galaxies peak at a redshift z ~ 15, while galaxies made of 
P opulation III stars peak a t a redshift z 11. This trend is in agreement with that shown in Fig. 1 
of lTrenti Stiavellil (|2009l ) once one considers that Pop II formation in that figure was for all halos 
regardless of their mass while here it is limited to halos of the prescribed virial temperature. The 
presence of a peak at z ~ 15 in the enriched galaxies is due to this mass constraint. 



3.2. Clustering Properties 

In Fig. [1] we show an analytical estimate of the bias as a function of redshift for halos that can 
contain Population III stars. Clearly the bias decreases rapidly with redshift as the Population III 
halos become more common. However, Fig. Q] does not distinguish between halos that are polluted 
by winds and those that are truly metal free (self-enrichment does not influence the bias). To address 
the expected clustering properties of the first galaxies we need to resort to the numerical simulation 
results, which of course carry spatial information on the presence of neighbouring galaxies capable 
of polluting pockets of metal free stars. This effect is important only for redshift z < 15 because 
at higher re dshift the metal o utflows do not have enough time to propagate far from their host 



galaxies (see lTrenti et al.ll2009l ) 



The cosmological simulation allows us to measure the bias for both the chemically enriched 
"first" galaxies population and the Population III galaxies. To quantify their clustering we use the 
positions of galaxy-hosting dark matter halos to construct the three dimensional two-point corre- 
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lation function 0(r) in comoving space, denned as the excess number of pair counts at separation 
r over those from a uniform random distribution. We generate a random un iform distribution of 
40000 points within the simulation box and adopt the lLandv fc Szalavl (|1993l ) estimator: 



l + 0(r) 



DD(r) - 2DR(r) + RR(r) 
RR{r) : 



(1) 



where DD(r) are the halo-halo pair counts with distance bound in an annulus centered on r. 
Similarly DR{r) are the halo-random pair counts and RR(r) the random-random counts within 
the same annulus. 

In Fig. [3] we show for redshift z = 9.5 the correlation function as a function of radius for 
chemically enriched protogalaxies, that is for all dark matter halos with 10 4 K < T v i r < 2xl0 4 K that 
had past star formation bursts (0p O pii marked as "Pop II Gal" , solid red) and for galaxies containing 
Population III stars, that is dark matter halos with T v i r > 10 4 K that are chemically pristine with 
respect to both wind and self-enrichment (0p O piii, blue lineal We also show the correlation function 
for extremely metal poor protogalaxies, at or above the critical metallicity Z ~ 10~ 3,5 Z , enriched 
by winds but otherwise not self-enriched, again hosted in halos with 10 4 K < T vir < 2 x 10 4 K ( green 
line). In the same figure we also s how the ratio of 6(r)p np iii/6(r)p np n to qua ntify the relative bias 
of these two galaxy populations (jPorciani et al.lll999l ; ISheth &: Tormenl ll999). Clearly for radii in 
excess of a few tens of kpc the chemically enriched galaxies exhibit stronger clustering than the 
galaxies containing Population III stars. The low bias for Population III galaxies compared to their 
metal enriched counterparts extends up to z ~ 15 (see Fig. [4]). 



4. Discussion and Conclusions 

Fig. [3] predicts that galaxies containing Population III stars will be essentially uncorrelated 
on scales larger than ~ 150 comoving kpc. At redshift z = 9.5, 1 arcsec corresponds to roughly 
50 comoving kpc in the transverse direction. This implies that the typical arcmin scales of the 
imaging instruments on the James Webb Space Telescope correspond to a few comoving Mpc and 
on those scales we should not detect any significant clustering for galaxies with primordial chemical 
composition. In contrast, enriched protogalaxies could be very clustered on small scales, as they 
live in high density regions. Thus, in an ultradeep field with JWST we expect to detect individual 
unevolved galaxies at z ~ 10, but small clusters of chemically evolved objects. From Fig. [2] we see 
that the rest frame rate of galaxy formation per comoving volume at z ~ 10 is ~ 2 x 10 -8 Mpc -3 
yr _1 (estimating the halo formation rate from the simulations gives a similar number). Assuming 
that the burst of Population III stars is short lived, we expect that the galaxy will be visible as 
Population III galaxy only for a time of the order of the stellar lifetime (assumed to be 2 x 10 6 yrs). 



lr There are no chemically pristine halos with T v i r > 2 x 10 4 K in our si mulation because we assume Population 
III star formation once halos are capable of Lya cooling ( Trenti et al. 201Clh 
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Considering that the comoving volume per unit redshift at z ~ 10 is ~ 2000 comoving Mpc 3 over 
an area of one arcmin squared, we predict a surface density of ~ 80 galaxies per square arcmin 
and per unit redshift. However, there are large uncertainties as to the luminosity such a galaxy 
would have. In a deep exposure, the James Webb Space Telescope should be able to detect a young 
cluster with a stellar mass of 10 6 mJ. Most of the pristine galaxies will have lower stellar masses. 
At a virial temperature T = 10 4 K at z ~ 10 a galaxy has a halo mass of ~ 3.7 x 10 7 M Q and only 
~6x 10 6 M in baryons. It is at this stage unclear whether the Population III stellar burst in a 
proto-galaxy would produce many stars and with what efficiency e, but e > 0.15 is unlikely in the 
first burst. Even if the majority of these objects will be too faint, their large number density and 
the expected lack of strong correlations are encouraging as they make it more likely that we could 
detect some of these objects through a lensing cluster. A Population III galaxy of 10 6 M of stars 
formed over a stellar lifetime forms s tars at about 0.3 M Q yr _1 and would have a Lyct luminosity 



of ~ 3 x 10 42 erg s 1 ( Schaererll2003l ) and a Lya line flux at z = 10 of ~ 1.9 x 10 18 erg s 1 cm 



These fluxes are close to what can be detected today in a len sing cluster after modest amplification 



and is in fact in agreement with some reports of detections (jStark et al.1 120071 ) 



The (chronologically) first galaxies at the same redshift and the extremely poor galaxies are 
instead expected to be highly correlated at least on scales of a few arcsec. Thus, to first order 
our models predict that one should observe small clusters of chemically evolved galaxies but only 
isolated primordial galaxies. 
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support from the University of Colorado Astrophysical Theory Program through grants from NASA 
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Fig. 1. — Left panel. Minimum halo mass for star formation (solid black line) required after 
accounting for H<2 photo-dissociation due to photons in the Lyman-Werner bands (figure from 



Trenti et al.l 120091 ) . At very high-redshift (z > 30) the gas can cool in small mini-halos, with 



T v i r ~ 10 3 K. As the global star formation rate increases, H2 cooling is suppressed and the minimum 
mass of a mini-halo capable of forming stars rises steadily. At z < 15 only halos with T v i r > 10 
K can form stars. Right panel: bias for halos at the minimum mass required for s tar fo rmation as 
shown in the left panel. The bias has been computed using the ISheth &: Tor men! ()1999l ) model. 
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Fig. 2. — Formation rate of proto-galaxies per comoving volume (in halos with 10 4 K < T v i r < 
2 • 10 4 K that are either chemically pristine (solid blue line) or with a low metallicity of the order of 
the critical one (red dashed line). The rate (re st-frame and comoving) has been obtained from an 
analytical Extended Press-Schechter model (see lTrenti Stiavellill2009l ) considering self-enrichment 
only, thus galactic metal outflows are neglected . If these were to be considered they would reduce 
the formation rate of PopIII proto-galaxies (see lTrenti et al.l l2009). The first galaxies to be formed 
at very high-redshift are very likely to be made of Population II (metal enriched) stars. 
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Fig. 3. — Left panel. Three-dimensional two-point correlation function for metal-free (blue) and 
metal-enriched (red) proto-galaxies, residing in halos with 10 4 K < T v i r < 2 • 10 4 K, at z = 9.5. We 
also show the correlation for extremely metal poor galaxies (green line) that have b een enriched 



by win ds. The figure has been obtained from the cosmological simulation presented in lTrenti et al 



(J2009j) and includes metal enrichment from outflows propagating at 60 km s . Lines are an average 
over 9 simulation snapshots with 9 < z < 10. Right panel: ratio of the Pop III to Pop II correlation 
functions, which further highlights that metal- free proto-galaxies are significantly less clustered on 
small scales (r < 1 Mpc/h) than their metal-enriched counterparts with similar halo mass. 
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Fig. 4. — Left panel. Three-dimensional two-point correlation function for metal-free (blue) and 
metal-enriched (red) proto-galaxies as in Fig. [3] measured at z = 14.5. The figure has been obtained 
by averaging over 5 simulation snapshots with 14 < z < 15. Right panel: ratio of the two-correlation 
functions. 



